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TECHNICAL  REPORT  SUMMARY 


Codas  of  local  earthquakes  from  Monticello  Reservoi r,  South 
Carolina,  and  Mammoth  Lakes,  California,  have  been  examined  to  determine 
scattering  parameters  and  attenuation  for  the  crust  in  these  areas.  The 
backscattering  turbidity  (or  backscattering  cross-section  per  unit 
volume)  and  the  total  coda  Q  have  been  determined.  In  both  areas, 
strong  scattering  is  observed  in  the  near  surface  crust  (upper  10  km), 
sufficiently  strong  that  multiple  scattering  is  occurring.  Because 
multiple  scattering  is  occurring,  coda  Q  at  short  times  will  under¬ 
estimate  Q  experienced  by  a  short  pulse.  At  long  times,  the  coda  at 
Monticello  indicates  the  presence  of  a  low-scattering  crustal  channel 
within  which  energy  probably  spreads  horizontally;  this  channel  is 
either  absent  or  more  strongly  scattering  at  Mammoth  Lakes.  Results  at 
both  sites  indicate  that  both  in  the  near  surface  region  and  in  the 
aforementioned  crustal  channel,  scattering  is  probably  an  important 
contributor  to  Q  (at  least  half  of  the  attenuation).  The  scatterers  are 
probably  predominantly  forward  scattering  (such  as  velocity  fluctua¬ 
tions),  and  in  the  near  surface  there  are  many  scatterers  of  size  50  m 
and  greater. 
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LONG  AND  SHORT  CODAS  IN  CALIFORNIA  AND  SOUTH  CAROLINA 

Introduction 

This  report  is  a  follow-up  report  to  an  earlier  study  (Dainty  and 
Tie,  1984)  in  which  the  local  coda  of  some  aftershocks  of  the  May  25  and 
27,  1980,  earthquake  near  Mammoth  Lakes,  California,  were  analyzed. 

This  analysis  indicated  that  for  coda  lengths  of  10  seconds  or  less,  the 
scattering  implied  by  the  decay  of  the  coda  was  not  consistent  with  the 
scattering  implied  by  the  amplitude  of  the  coda  relative  to  the  ampli¬ 
tude  of  the  direct  S  wave.  Briefly,  the  scattering  implied  by  the  coda 
amplitude  should  lead  to  more  rapid  decay  of  the  coda  than  is  observed 
at  frequencies  around  3  Hz,  if  the  single  scattering  theory  used  to 
analyze  the  data  is  correct. 

I  have  examined  other  local  earthquake  records  to  see  if  the 
results  presented  above  are  correct  for  these  codas.  Besides  records  of 
length  10  seconds  or  less,  I  have  analyzed  records  of  substantially 
longer  length  to  see  if  the  effects  are  seen  at  longer  times.  Different 
records  from  Mammoth  Lakes  (Archuleta  et  al_. ,  1982),  analyzed  by  a 
different  investigator,  and  records  from  Monticello  Reservoir,  South 
Carolina  (Fletcher,  1982),  will  be  used.  These  are  the  records  analyzed 
by  Duckworth  (1983)  and  Dainty  and  Duckworth  (1983). 

Theoretical  Summary:  Single  Scattering  and  Multiple  Scattering 

The  standard  theory  used  to  analyze  codas  is  the  single  scattering 
theory  due  to  Aki,  summarized  in  Aki  and  Chouet  (1975).  This  theory  has 
been  used  in  earlier  reports.  The  theory  assumes  a  whole  space  in  which 
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scatterers  are  distributed  randomly.  These  scatterers  may  be  approxi¬ 
mated  by  some  body  of  specified  geometrical  shape  or  may  be  defined  by 
their  autocorrelation  functions.  Since  an  elastic  medium  has  (at  least) 
three  independent  relevant  material  constants,  two  elastic  constants  and 
density,  the  three  quantities  must  be  specified,  although  we  expect  some 
correlation  between  the  three  quantities.  Generally,  it  has  been 
assumed  that  only  S  to  S  scattering  is  important,  but  recent  results  by 
Wu  and  Aki  (1984)  suggest  that  S  to  P  scattering  could  also  be  impor¬ 
tant.  P  to  P  and  P  to  S  scattering  would  be  less  important  because 
earthquakes  do  not  generate  as  much  P  as  S. 

Briefly,  the  theory  alluded  to  implies  that  at  long  times  the  coda 
power  spectrum  shall  decay  with  time  t  after  origin  as 


Pc(“.t)  =  Ab(u)t  exp(-wt/Q) 

AbU)  *  P$U)  .  exp(u)t$/Q)  .  8wV  .  g(ir,oj)  .  t‘ 


1/Q  -  1/C^  +  1/QS 


1/Q-  =  G(u)  V/t 


P  (ui)  is  the  square  of  the  S  wave  spectrum,  V  is  the  seismic  velocity 
(assumed  to  be  3.2  km/sec,  typical  of  S  waves  in  the  crust,  in  this 
study),  t  is  the  S  travel  time,  Q..  is  the  Q  due  to  anelasticity ,  and  Q$ 
is  the  Q  due  to  scattering.  The  scattering  is  described  by  two  para¬ 
meters.  The  backscatteri ng  turbidity  g(ir,u)  or  backscattering  cross 
section  per  unit  volume  describes  the  strength  of  backscattering  and 
controls  the  amplitude  of  the  coda  relative  to  the  amplitude  of  the  S 
wave,  as  seen  in  (2).  The  total  turbidity  G(w),  or  total  cross-section 
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per  unit  volume,  describes  the  strength  of  total  scattering  in  all 
directions  and  controls  the  Q  due  to  scattering  through  (4).  Another 
parameter  that  I  have  calculated  is  the  apparent  total  turbidity 


Ga(.)  .ft 


If  the  intrinsic  attenuation  can  be  neglected,  the  apparent  total 
turbidity  is  the  same  as  the  total  turbidity.  Generally,  in  single 
scattering. 


G,(w)  >.  G(,J))  >.  g(n»al) 


So  far,  the  theoretical  discussion  has  assumed  that  the  coda  con¬ 
sists  of  scattered  body  waves  of  one  type  in  a  whole  space.  In  view  of 
the  results  to  be  presented  casting  doubt  on  this,  some  other  possibili 
ties  must  be  discussed.  The  simplest  is  the  case  of  singly  scattered 
surface  waves,  rather  than  body  waves,  from  a  near  surface  source: 
Equation  (1)  becomes 


Pc(oj,t)  =  As(w)t“  exp(-wt/Q) 


ASU) 


P;(w)  .  exp(wt^/Q)  .  2irV  .  g' (*,<*>) 


1/Q  =  1/Q!  +  G' (w)V'/t 


The  primed  quantities  in  (8)  and  (9)  indicate  that  these  variables  have 
the  same  definitions  as  before,  but  for  surface  waves  rather  than  S 
waves.  Another  possibility  is  of  two  or  more  singly  scattered 
"channels"  for  the  transmission  of  energy  that  appears  in  the  coda--for 
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example,  one  channel  could  be  body  wave  scattering  and  the  other  surface 
wave  scattering.  Then  the  coda  power  spectrum  would  become 

PcKt)  =  PxKt)  +  P2(u.,t)  (10) 

where  and  P2  are  the  contributions  from  each  channel. 

Another  important  issue  is  that  of  multiple  scattering  as  opposed 
to  single  scattering.  Dainty  and  Toksoz  (1977)  demonstrated  that  in 
the  limit  of  extreme  multiple  scattering,  seismic  energy  diffuses.  The 
coda  power  spectrum  for  a  -point,  impulsive  source  in  a  whole  space  then 
becomes  (Dainty  and  Toksoz,  1981): 

PcUt)  =  Eb(u)t~3/2  exp[-3R2GH/(4Vt)]  exp[-o)t/Qi ]  (11) 

— >  Eb(u)t‘3/2  e[-wt/Qi ]  ,  R  0  (12) 

Eb(.)  =  Ps(oj)  .  /27G3(.)  /  (8 J7?)  (13) 

Here  R  is  the  source-receiver  distance,  taken  to  be  small  in  (12).  Note 
that  (12)  is  rather  similar  to  (1)  or  (7)  in  that  it  indicates  that  at  a 
fixed  frequency,  for  short  source-receiver  distance  (the  usual  experi¬ 
mental  situation),  the  coda  power  spectrum  should  decay  with  time  as  trie 
product  of  an  inverse  power  law  intermediate  between  (1)  and  (7)  and  an 
exponential.  Thus  if  an  attempt  was  made  to  fit  equation  (1)  or  (7)  to 
a  situation  in  fact  described  by  equation  (12),  an  apparently  good  fit 
would  be  obtained.  However,  the  interpretation  of  the  parameters  would 
change.  In  (12),  the  exponential  decay  gives  the  value  of  ,  the 
intrinsic  attenuation,  and  not  Q  as  defined  by  (3)  or  (9).  In  fact, 

Qi  >  Q 


(14) 


I  will  assume  that  the  general  effect  of  multiple  scattering  is  to 
increase  the  value  of  Q  measured  as  compared  to  that  implied  by  (3)  or 
(9).  This  will  decrease  G  (^)  calculated  from  (5). 

d 

In  this  paper  g(w,w)  has  been  calculated  from  (2)  or  (8),  i.e., 
under  the  assumption  of  single  scattering.  Essentially,  this  means  that 
gU.jj)  has  been  calculated  from  the  ratio  of  the  coda  power  to  the  S 
wave  power,  with  a  (variable)  factor.  From  (13),  this  ratio  should  be 
large  for  the  case  of  multiple  ("strong")  scattering,  since  G( 1  '  *11  be 

large;  thus  large  values  of  g(ir,u)  should  be  found.  Combir  ,  this 
result  with  the  result  for  G  (w)  just  discussed,  increasing  g  -  and 
simultaneously  decreasing  G  (m)  should  be  indicative  of  the  onset  of 

a 

multiple  scattering. 

Other  criteria  for  determining  if  multiple  scattering  is  occurring 
have  been  proposed.  Dainty  and  Toksoz  (1981)  give 

1/Q.  «  GV/u>  (15) 

as  a  sufficient  but  not  a  necessary  condition  for  multiple  scattering. 
Sato  (1977),  on  the  other  hand,  gives 

tm  >  1/(GV)  (16) 

as  a  necessary  but  not  a  sufficient  condition  for  multiple  scattering, 
where  t  is  the  total  length  of  the  coda  after  origin. 

The  remaining  theoretical  question  is  the  characterization  of  the 
scatterers,  which  controls  the  relationship  between  g(ir,u>)  and  G ( J] . 
Previously,  it  had  been  hoped  that  the  measurements  reported  here  would 
answer  this  question.  However,  whilst  they  appear  to  shed  some  light  on 


the  subject,  it  appears  that  there  are  too  many  parameters  involved  to 
allow  unique  i nterpretati ons.  Briefly,  two  types  of  models  have  been 
used  to  characterize  the  scatterers--either  bodies  of  a  specified  shape, 
distributed  randomly  (e.g..  Dainty,  1981),  or  random  fluctuations  of  the 
medium  described  by  their  autocorrelation  (e.g.,  Wu,  1982).  There  are 
some  similarities  between  the  various  models--if  there  is  a  length  scale 
a  associated  with  the  scatterers  (the  object  size  if  the  scatterers  are 
modelled  as  bodies,  or  a  correlation  scale  length  for  random  fluctua¬ 
tions),  then  for  wavelengths  A  much  longer  than  this  length  scale  the 
scattering  is  weak  but  strongly  frequency  dependent,  increasing  as  the 
fourth  power  of  frequency  (Rayleigh  scattering).  The  scattered  inten¬ 
sity  increases  until  a/A  ^  1.  For  media  with  sharp  boundaries,  back- 
scattering  is  approximately  frequency  independent  for  a/A  >>  1  (geo¬ 
metrical  scattering;  Dainty,  1984),  or  declines  rapidly  if  there  are  no 
sharp  boundaries.  Thus,  if  a  variety  of  length  scales  of  scatterers  are 
present  in  the  earth,  the  length  scales  that  are  the  same  order  of 
magnitude  of  the  wavelength  or  greater  are  likely  to  be  important. 

An  important  area  of  difference  between  the  models,  however,  lies 
in  the  predicted  dependence  of  the  scattered  intensity  on  the  angle 
between  the  original  wave  and  the  scattered  wave.  Almost  all  possi¬ 
bilities  seem  to  be  allowed,  including  isotropic  scattering  (spheres  at 
high  frequency),  strong  forward  scattering  (acoustic  velocity  fluc¬ 
tuations  at  high  frequency;  Dainty,  1984),  or  strong  backscatteri ng 
(impedance  fluctuations;  Wu  and  Aki,  1984).  At  any  given  frequency,  the 
angular  dependence  of  scattering  controls  the  ratio  of  the  back- 


8 


scattering  turbidity  g(ir,uO  to  the  total  turbidity  G(uj).  For  the  case 
of  isotropic  scattering 

gU,a))/G(u)  =  1  / ( 4ir )  10%  (17) 

Ratios  of  g(ir,uj)  to  G(jj)  substantially  less  than  0.1  (10%)  would  indi¬ 
cate  a  medium  that  predominantly  scatters  in  the  forward  direction; 
Dainty  (1984),  for  an  acoustic  medium  with  velocity  fluctuations  at  high 
frequencies,  obtains 


g(ir,u>)/G(u>)  =  l/(56ir)  n,  1%  (18) 

Likewise,  ratios  of  g(ir,co)  to  G(oj)  substantially  greater  than  0.1  would 
indicate  a  predominantly  backscattering  medium.  Wu  and  Aki  (1984) 
suggest  that  an  elastic  medium  with  impedance  fluctuations  would  have 
these  properties. 

Analysis  and  Results 

The  data  used  in  this  paper  are  mainly  those  of  Duckworth  (1983), 
with  data  from  Dainty  and  Tie  (1984)  used  for  comparison  purposes. 
Analysis  methods  are  given  in  Duckworth  (1983)  and  Dainty  and  Duckworth 
(1983).  Briefly,  codas  are  analyzed  by  determining  the  power  spectrum 
as  a  function  of  time  in  a  moving  window  1.28  sec  long  for  the  results 
of  Duckworth  (1983),  or  0.64  sec  for  the  results  of  Dainty  and  Tie 
(1984).  The  power  spectra  are  averaged  over  octave  bands  and  then 
fitted  as  a  function  of  time  at  a  specified  center  frequency  to  equation 

(1)  or  equation  (7)  to  obtain  Q  and  A^  or  A  .  Then  equations  (5)  and 

(2)  or  (8)  as  appropriate  were  used  to  calculate  G  (j)  and  g(ir,u>). 

a 

Further,  following  Duckworth  (1983),  the  coda  was  split  into  two  time 
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sections:  times  less  than  10  seconds  after  origin  and  times  greater 
than  10  seconds.  This  followed  visual  examination  of  plots  of  the 
Monticello,  South  Carolina,  codas  as  power  at  a  frequency  as  a  function 
of  time.  There  appeared  to  be  a  change  in  the  decay  rate  of  the  coda  at 
about  10  seconds  after  origin,  and  separate  fits  were  made  to  examine 
this,  both  on  the  Monticello,  South  Carolina,  records  and  the  Mammoth 
Lakes,  California,  records. 

The  records  used  are  digital  records  used  by  Fletcher  (1982)  at 
Monticello  Reservoir,  South  Carolina,  and  by  Archuleta  et  aQ_.  (1982)  at 
Mammoth  Lakes.  The  digital  recordings  were  kindly  supplied  by  Paul 
Spudich  of  the  United  States  Geological  Survey.  The  seismograms  are 
sampled  at  a  rate  of  200  samples  per  second  with  a  12  bit  sample.  The 
antialiasing  filter  has  a  corner  of  50  Hz;  analysis  has  been  limited  to 
frequencies  less  than  this.  The  events  and  the  stations  used  in  this 
study  are  shown  in  Figures  1  and  2.  Table  1  indicates  which  events  and 
stations  were  used  in  such  analysis. 

Results  of  the  analysis  are  shown  in  Figures  3  through  6.  For  each 
seismogram  analyzed,  both  equation  (1)  and  equation  (7)  were  fit.  The 
results  from  equation  (1)  will  be  called  the  "spherical"  case  and  the 
results  from  equation  (7)  the  "cylindrical"  case.  Figures  3  through  6 
show  average  values  for  both  G  (.»)  and  g ( -rr , jo )  for  each  fit  at  Mammoth 

d 

and  Monticello  for  long  (>10  sec)  and  short  (<10  sec)  codas.  Note  that 
results  from  the  "spherical"  case  usually  give  lower  numerical  values 
than  results  from  the  cylindrical  case,  but  the  general  trend  of  the 
data,  and  the  relative  values  of  G  (.a)  and  g(n,w)  are  similar.  Earlier 

d 

results  of  Dainty  and  Tie  (1984)  using  different  events  and  a  shorter 
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I  able  lb.  Locations  of  earthquakes  and  recording  stations  at 
Mammoth  Lakes,  California. 
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Figure  1.  Recording  stations  and  epicenters  at  Monticello  Reservoir 
South  Carolina.  Stations  -  1,  epicenters  -  * 


Recording  stations  and  epicenters,  Mammoth  Lakes,  California 
Stations  -  A,  epicenters  -  *.  Outline  of  Long  Valley 
Caldera  and  major  faults  shown. 


,  closed  symbols,  cylindrical  spreading  assumed.  Circles 


turbidity  from  Dainty  and  Tie  (1984);  open  hexagons,  backscattering  turbidity.  Dainty 
and  Tie  (1984). 


window  are  shown  in  Figure  5  for  comparison;  these  earlier  results  refer 

to  a  spherical  case,  short  coda,  and  were  obtained  from  Mammoth 

seismograms.  These  results  will  be  discussed  in  the  next  section. 

Sometimes,  negative  values  of  Q  were  obtained  when  fitting  equation 

(1),  and  less  frequently  equation  (7).  This  indicates  that  the  coda 

_2 

power  at  that  frequency  is  decaying  less  rapidly  than  t  if  equation 
(1)  is  being  fit,  or  less  rapidly  than  t~*  if  equation  (7)  is  being  fit. 
Such  values  have  been  counted  as  Q  »  (G  (w)  -4  0)  in  the  averages. 

a 

In  the  log-log  plots,  values  of  G  (w)  correspondi ng  to  this  case  in  the 

a 

average  have  been  plotted  at  the  lowest  scale  values. 

Discussion 

The  first  question  that  will  be  addressed  is  the  issue  of  spherical 
(body  wave)  spreading  as  opposed  to  cylindrical  (surface  wave)  spreading 
of  the  scattered  waves  in  the  coda.  In  Figures  3  through  6,  the  results 
of  spherical  spreading  fits  are  shown  as  open  symbols  and  the  results  of 
cylindrical  spreading  fits  as  closed  symbols.  It  appears  that  in  most 
cases,  an  equally  good  fit  can  be  obtained  using  either  assumption.  The 
major  exception  to  this  statement  is  seen  in  Figure  4,  Monticello  codas 
more  than  10  seconds  long,  where  negative  values  of  Q  are  obtained  for 
frequencies  near  10  Hz  for  spherical  spreading  fits.  This  suggests  that 
cylindrical  spreading  along  a  "channel"  is  responsible  for  the  coda  at 
Monticello  after  10  seconds.  Because  the  effect  is  noticeable  at  times 
as  short  as  10  seconds,  the  channel  must  exist  in  the  crust  within  about 
15  km  of  the  surface,  since  the  scattered  energy  must  get  from  a  shallow 
(%1  km  deep)  source  to  the  scatterers  and  back  to  the  surface  receiver 
in  as  little  as  10  seconds. 


It  is  clear  from  the  results  obtained  and  shown  in  Figures  3 

through  6  that  multiple  scattering  is  occurring  at  frequencies  less  than 

about  20  Hz  in  all  cases  except  the  aforementioned  Monticello  codas  at 

times  longer  than  10  seconds.  For  all  these  cases  there  is  a  rise  of  the 

backscattering  turbidity,  gU.w)  for  frequencies  less  than  20  Hz, 

without  a  corresponding  rise  in  G  (w).  Indeed,  in  the  case  of  Mammoth 

Lakes  (Figures  5  and  6),  g(iT,aj)  is  greater  than  G  (a)  for  the  lowest 

8 

frequency  studied,  3  Hz.  This  violates  equation  (6)  severely,  and  is 
not  possible  under  the  assumption  of  single  scattering.  A  similar 
result  was  obtained  by  Andrews  (1982).  Furthermore,  in  the  case  of 
Mammoth  Lakes  the  apparent  total  turbidity  declines  by  at  least  a  factor 
of  two  at  3  Hz,  suggesting  that  at  least  half  of  the  decay  at  higher 
frequencies  is  due  to  scattering.  The  case  of  Monticello  codas  at  times 
longer  than  10  seconds,  however,  seems  to  fit  the  single  scattering 
theory  well,  obeying  condition  (16). 

The  nature  of  the  scatterers  may  be  partly  determined  by  the  data 
presented  here.  Basically,  the  approach  is  to  compare  the  relative 
values  of  g(ir,w)  and  G  (u>).  Consider  first  Figure  4,  the  data  from 
Monticello  for  codas  longer  than  10  seconds.  The  backscattering  tur¬ 
bidity  g(u,w)  and  the  apparent  total  turbidity  G  (w)  show  an  excellent 
correlation  with  each  other.  When  the  apparent  total  turbidity 
decreases  around  10  Hz,  so  does  the  backscattering  turbidity.  This 
suggests  that  the  apparent  total  turbidity  is  mainly  due  to  scattering 
and  thus  a  close  approximation  to  the  true  total  turbidity  G(a>).  The 
value  of  the  backscattering  turbidity  is  about  1%,  or  less,  than  the 
total  turbidity,  which  indicates  a  medium  that  mainly  forward  and/or 


20 


side  scatters,  such  as  a  velocity  fluctuation  (Dainty,  1984;  Wu  and  Aki, 
1984). 

The  other  cases  examined  (Figures  3,  5,  and  6)  are  not  so  easy  to 
interpret  for  two  reasons.  One  is  the  effect  of  multiple  scattering 
already  noted,  preventing  a  measure  of  g ( it , -^ )  and  G(w).  This  is  the 
most  serious  difficulty  at  low  frequency.  At  high  frequencies,  we  note 
from  Figure  7  that  the  back  scattering  turbidity  g(ir,o)  has  the  same 
behavior  for  all  three  cases,  quite  distinct  from  the  Monticello  codas 
longer  than  10  seconds.  From  Figure  7,  however,  the  apparent  total 
turbidity,  G  (u),  at  Mammoth  Lakes  (both  cases)  seems  to  drop  at 
frequencies  above  20  Hz.  This  may  be  due  to  contamination  by  noise, 
which  would  tend  to  make  the  coda  decay  seem  less  than  it  really  is, 
without  affecting  the  determination  of  the  backscatteri ng  turbidity 
g(it,'j),  as  suggested  by  Figures  7  and  8.  In  support  of  this,  note 
(Andrews,  1984,  personal  communication)  that  there  is  very  little  energy 
in  the  Mammoth  Lakes  sources  above  20  Hz. 

If  in  tact  the  cases  of  the  Monticello  codas  less  than  10  seconds 
long  and  the  Mammoth  Lakes  codas  (both  coda  lengths)  are  similar,  as 
suggested  by  Figure  8,  we  may  compare  the  two  Monticello  cases  of  codas 
less  than  10  seconds  long  (Figure  3)  and  greater  than  10  seconds  long 
(Figure  4).  Figures  7  and  8  may  also  be  used.  We  see  that  both  g(ir,w) 
and  G  (.jj)  are  a  factor  of  5  to  10  larger  for  the  short  codas  as  compared 

a 

to  long  codas  for  frequencies  greater  than  20  Hz.  This  suggests  that 
the  channel  in  which  short  coda  energy  propagates  is  similar  to  that  for 
long  codas  in  its  scattering  properties,  but  scatters  more  strongly.  At 
frequencies  below  20  Hz,  little  can  be  said.  It  is  possible,  through 


equation  (15),  for  the  scattering  properties  at  low  frequencies  to  be 
the  same  as  those  at  high  frequencies,  and  yet  for  there  to  be  multiple 
scattering  at  low  frequencies  (<20  Hz)  and  not  at  high  frequencies  (>20 
Hz),  if  is  around  1000.  However,  this  is  the  approximate  value  of 
the  total  Q  at  high  frequency  (>20  Hz),  leaving  no  room  for  a  scattering 
contribution.  Accordingly,  it  seems  more  likely  that  there  is  a  group 
of  scatterers  of  size  50  m  and  larger  causing  the  increased  scattering 
implied  by  the  rise  in  g(ir,u)  seen  for  codas  of  less  than  10  seconds 
length  at  Monticello. 

Finally,  what  is  the  nature  of  the  difference  in  coda  observed 
before  and  after  10  seconds  at  Monticello,  South  Carolina?  I  have 
interpreted  this  in  terms  of  equation  (10),  i.e. ,  as  two  different 
regions  of  the  crust  which  scatter  differently.  Note  that  if  one 
channel  has  high  scattering,  it  will  tend  to  dominate  at  short  times 
because  g(ir,co),  controlling  the  excitation,  is  large,  but  will  be 
negligible  at  long  times  because  G(a)),  controlling  the  decay,  is  large. 
Thus  the  short  time  (<10  sec)  coda  shows  the  stronger  scattering.  Since 
this  coda  occurs  immediately  after  the  S  wave,  and  at  Mammoth  Lakes  does 
not  appear  to  be  dependent  on  focal  depth  (Tie,  personal  communication), 
I  suspect  that  it  is  near  surface,  spherical  (?)  scattering.  The  long 
coda  (>10  sec)  may  represent  cylindrical  (?)  propagation  in  a  crustal 
"channel",  possibly  at  somewhat  deeper  depths  (10  to  15  km?). 

From  Figure  8  especially,  it  seems  that  there  is  little  difference 
between  the  codas  at  Mammoth  Lakes,  both  long  and  short,  and  the  short 
(<10  sec)  coda  at  Monticello.  This  indicates  that  the  upper  crust  is 
rather  similar  at  these  two  tectonically  very  different  sites.  However, 


at  times  greater  than  10  seconds  the  Monticello  coda  shows  evidence  of  much 
weaker  scattering  than  at  times  less  than  10  seconds,  whereas  the 
Mammoth  Lakes  data  show  little  if  any  significant  difference  between 
codas  in  the  two  time  ranges  (Figure  8).  Either  the  crustal  low- 
scattering  channel  seen  at  Monticello  is  not  present  at  Mammoth  Lakes 
(the  most  likely  conclusion)  or  it  scatters  as  strongly  as  the  near¬ 
surface  channel.  Because  propagation  in  the  crustal  low-scattering 
channel  is  suspected  to  be  "cylindrical",  i.e. ,  confined  to  the  hori¬ 
zontal  direction,  it  may  be  the  same  channel  as  that  used  by  Lg.  If 
this  is  true,  the  results  obtained  have  predicted  that  Lg  would  propa¬ 
gate  more  efficiently  at  Monticello  than  at  Mammoth  Lakes,  a  prediction 
that  is  certainly  true  for  the  wider  regions  in  which  these  areas  are 
located  (Nuttli,  1973).  Further,  the  Q  measured  from  long  time  codas 
would  be  representati ve  of  Lg  Q  (Mitchell  and  Nuttli,  1983). 
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